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The  Scalar  Riemann  Problem  in  Two  Spatial  Dimensions: 
Piecewise  Smoothness  of  Solutions  and  its  Brealtdown 

W.  B.  Lindqidst  ^ 

Courant  Institute  of  Mathonatical  Sciences 

New  York  University 

New  York,  N.Y.  10012 


ABSTRACT 


Consider  the  scalar  quasilinear  equation 

duU,x)    .    «  a//(w(f,  x))  _ 
dt       ^i,       dx, 

for  n  =  1,  2  ,  f,  d  C^  -.R  "R.  For  n  =  2,  we  define  the  two  dimensional  Riemann 
problem  and  show  the  unique  (in  the  sense  of  Knizkov)  solutions  are  piecewise  smooth  for 
/i  ■  /j  ■  /,  /  purely  convex  or  having  a  single  inflection  point.  A  mechanism  leading  to  a 
presumed  loss  of  piecewise  smoothness  is  presented  for  /  having  three  ot  more  inflection 
points.  The  analysis  is  based  on  a  study  of  the  generalization  of  the  one  dimensional 
Riemann  problem  to  allow  for  initial  data  having  a  finite  number  of  jump  discontinuities 
with  constant  data  or  rarefaction  waves  between  jumps. 


*  Supported  by  the  Applied  Mathematical  Sciences  subprogram  of  the  Office  of  Energy  Research, 
U  S  Dept.  of  Energy,  Contraa  DE-Aa)2-76ERO30T7. 


-2- 

1.  IntrododfaHi 

Consider  the  Caucfay  problon  for  the  scalar  quasilinear  equation 


SuiLA  ^  ^SfMilJlll  =  0  ,       n  =  l,2, 


at        (.1      dx, 

«(0,x)  =  «oW.  (11) 

with  f,:R-'R  ^  C^.  For  two  spatial  dimensions,  we  define  the  Rionann  problem  and  show  the 
unique  (in  the  sense  of  Kruzkov)  solutions  are  pieoewise  smooth  for  /i  "  /j  *  />  /  having  at  most 
one  inflection  point.  The  assumption  f^  =  f^  appears  to  be  natural  from  the  point  of  view  of 
ai^lications.  This  sufficiency  condition  follows  from  a  detailed  analysis  of  a  generalization  of  the 
Riemann  problem  in  one  dimension  to  allow  for  initial  data  having  a~  finite  number  of  jump 
discontintuities  where  the  data  between  points  of  jump  is  allowed  to  contain  rarefaction  waves  as 
defined  below.  For  these  generalized  one  dimensional  Riemann  problems,  we  show  the  unique  (in 
the  sense  of  Kruzkov)  solutions  are  piecewise  smooth  if  the  initial  data  has  a  single  jump 
discontinuity  and  /^  is  restricted  to  at  most  two  inflection  points.  If  the  initial  data  has  more  than 
a  single  jump  discontinuity  we  show  that  restricting  /^  to  at  most  one  inflection  point  guarantees  a 
piecewise  smooth  solution. 

We  present  a  construction  based  on  a  function  /^  having  three  inflection  points  which  shows 
a  possible  breakdown  of  pieoewise  smoothness  for  the  associated  one  dimensional  Riemann 
problem.  The  implication  of  this  example  is  that  the  scalar  Rionann  problem  in  two  dimensions 
can  presumably  fail  to  have  pieoewise  smooth  solutions  for  non  convex  /2  ■  /j  having  three 
inflection  points. 

A   bounded    measurable   function    u(r,  jc)    is    a    weak   solution    of    (1.1)    in   the   strip 

Uj.  -  [o,r]  X  /?"  if 


J  Jdtdx 


dt        <.i  dx, 


+      JdxUoix)^iO,x)     =0,  (1.2) 


1-0 


for  all  «})  €  Co*  (Hj-).  A  function  is  a  weak  solution  in  the  large  if  it  is  a  weak  solution  for  all  T. 
For  initial  data  Uq(x)  in  the  class  of  bounded,  measurable  functions,  Kruzkov  [2]  has  shown  the 
existence  and  uniqueness  of  a  boimded  measurable  function  which  is  a  weak  solution  to  (1.1)^ 


'  Kruskov's  proof  hdds  for  n  a  1  spatial  dimensions.  Here  only  n  s  2  corcems  us. 


The  uniqueness  (entropy)  condition  satisfied  by  this  weak  solution  is  [2] 

fjdtdx  sign(«-Jfc)|  (u  -  k)^  +  i  (/;(«)  -  /,(*))|^  1^0,  (1.3) 

for  any  real  constant  k  and  any  <(>  C  Cq'CUj-)  such  that  <)>  ^  0.    This  uniqueness  condition 
characterizes  the  allowed  discontinuities  in  the  weak  solution. 

Define  a  regular  point  of  the  weak  solution  u  as  a  point  p  for  which  du/dt  and  du/dx^  exist 
and  are  continuous  in  a  ndghbourhood  of  p,  and  at  which  the  differential  form  of  (1.1)  is 
satisfied.  The  function  u  shall  be  termed  smooth  in  an  open  neighbourhood  M  if  each  point  p  (.  M 
is  a  regular  point.  Let  a  smooth  surface  of  codimension  i  in  R"  x  R"^  be  a  surface  at  each  point  p 
of  which  a  unique  (modulo  minus  signs)  i-dimensional  space  of  normals  exist  and  for  which  the 
normal  space  varies  continuously  as  the  point  p  is  varied  over  the  surface.  A  piecewise  smooth 
surface  is  defmed  by  induction  on  the  codimension  i,  as  the  union  of  a  finite  number  of  smooth 
surfaces  with  piecewise  smooth  boimdaries  of  codimension  i+1;  each  of  the  finite  number  of 
pieces  being  everywhere  the  boundary  of  the  same  codimension  i  surfaces.  A  piecewise  smooth 
solution  u  €  R"  X  R*^  b  defined  as  a  function  that  is  smooth  except  on  a  discontinuity  set  that 
consists  of  piecewise  smooth  surfaces.  Apoint  cf  jump  is  a  point  p  of  discontinuity  of  a  function  u 
for  which  the  discontinuity  is  a  simple  jump  in  u  •  ie.  locally  the  point  ;>  lies  on  a  smooth 
codimension  1   surface  of  discontinuity  and  the  one  sided  limits  u~  =  lim  u(p  —  en)  and 

f-O* 

u"^  =  lim  u(p  +  en),  where  n  is  the  unit  normal  to  the  smooth  discontinuity  surface  at  p,  exist  at 

p.  A  singular  point  of  u  which  is  not  a  point  of  jump  shall  be  termed  an  irreguiar  point. 

Using  integration  by  parts,  the  arbitrariness  of  4*  and  appropriate  choices  of  the  constant  k, 
it  can  be  shown  that  (1.3)  is  equivalent  to  two  conditions  for  pants  of  jump  p  in  whose 
neighbourhood  the  unique  weak  solution  u  is  piecewise  smooth: 

n[«--«-,  /(U-) -/(«-)]  =  0,  (1.4) 

and 

n[k-u-,  m-f(u-)]^0,  (1.5) 

where  /  stands  for  the  vector  [/i,/2  ].    In  (1.5)  n  is  oriented  such  that  u~  s  u'*'  and  k  is  any 
constant  such  that  u"  s  Jk  s  u"^.   As  a  consequence  of  (1.4),  (1.5)  also  holds  with  u~  replaced  by 


u*.  Equation  (1.4)  will  be  denoted  the  jump  condition  for  the  discontinuity  at  p  and  (1.5)  the 
entropy  condition.  In  R^  x  R"^  (1.4)  is  the  familiar  Ranldne-Hugoniot  condition  and  (1.5)  is 
equivalent  to  the  entropy  condition  of  Qeinik  [5].  The  advantage  of  equations  (1.4)  and  (1.5) 
over  (1.3)  is  that  the  former  are  local  conditions  which  will  allow  construction  of  pieoewise  smooth 
global  solutions  by  piecing  together  local  solutions  which  individually  obey  (1.4)  and  (1.5). 

A  weak  solution  u  €  A"  x  /?'*'  to  (1.1)  which  is  piecewise  smooth  and  whose  discontinuity 
sets  are  pieoewise  smooth  surfaces  of  points  of  jump  obeying  (1.4)  and  (1.5)  then  satisfies  (1.3) 
and  is  therefore  the  unique  Kruzkov  solution  to  (1.1). 

The  existence  and  uniqueness  theory  gives  little  insight  into  the  actual  form  of  the  weak 
solution.  In  R^  x  R'^,  insight  is  provided  by  analysis  of  the  Riemann  problem  whose  solution  is 
piecewise  smooth  and  can  be  derived  in  terms  of  a  non-linear  wave  analysis.  It  can  be  shown  (see, 
for  example,  Lax  [3]  and  the  references  therein)  that  these  non-linear  waves  consist  of  rarefaction 
and  shock  waves.  The  initial  motivation  of  this  paper  was  to  consider  the  natural  generalization  of 
the  Riemann  problem  in  two  dimensions  and  investigate  its  solution  in  terms  of  two  dimensional 
waves  which,  by  implication,  is  a  restriction  to  the  set  of  piecewise  smooth  solutions.  In  this  paper 
we  determine  a  sufficient  condition  for  which  the  solution  to  the  two  dimensional  Riemann 
problon  is  piecewise  smooth.  In  a  companion  papa  [4]  we  investigate  the  form  of  these  solutions 
in  terms  of  non-linear  waves.  For  the  case  /i  =  /2  ■■  /.  the  sufficiency  condition  consists  of 
restricting  /  to  at  most  one  inflection  point.  We  conjecture  that  if  the  number  of  inflection  points 
of  /  is  at  least  three,  the  solution  to  the  generalized  Riemann  problon  defined  in  1  dimension  (and 
hence  for  the  Rionann  problon  in  two  dimensions)  may  fail  to  be  pieoewise  smooth.  Our 
argument  for  this  is  presented  in  52. 

The  generalizations  of  the  1 -dimensional  Riemann  problon  to  allow  multiple  jun^  and 
restricted  smooth  variation  in  the  initial  data  are  presented  in  S§2  and  3  .  In  S4,  the  two 
dimensional  Riemann  problem  is  defined  and  a  sufficiency  condition  for  piecewise  smooth 
solutions  obtained. 


2.  The  l-Dfanensioiial  Riemaiin  Problem 

Throughout  this  paper  the  notations  q'(p)  or  (q(p))'  shall  denote  d  q(p)  /  dp  for  any 
function  9  of  a  single  variable />. 

We  first  discuss  the  familiar  results  for  the  one  dimensional  Riemaim  problem  to  introduce 
notation  and  reformulate  results  in  a  manner  that  will  allow  generalization.  The  scalar  quasilinear 
eqiiation  (1.1)  in  one  spatial  dimension  (f^  «  /,  x^  *  x)  with  initial  data 

<^'-)  =  [    uUorx>xl:  (2.1) 

defines  the  1 -dimensional  Riemann  problem. 

The  Riemann  problem  (1.1),  (2.1)  is  invariant  under  the  similarity  transformation 
(f,  X  -  xq)  -  {ct,  c(x  -  xq)),  c  >  Q  ?  Hence  the  solution  is  constant  along  the  straight  lines 
(x  -  Xo)/r  =  constant.  In  a  neighbourhood  of  a  point  in  which  the  Riemann  problon  solution  is 
smooth,  (1.1)  implies  that  the  solution  will  be  constant  along  the  straight  lines 

=  /'(u),      t^t,,  (2.2) 


t-tc 


which  are  denoted  characteristics.  Characteristics  shall  be  denoted  («,  r^,  x^)  where  (r^,  x^)  is  the 
origin  point  of  the  characteristic  given  by  (2.2).  For  the  Rionann  problan  (1.1)  (2.1),  r^  =  0  f or 
all  characteristics. 

Let  p  =  (tp,  Xp)  be  a  point  on  a  characteristic  (m^,  fj,  Xj).  Let  n  be  a  normal  to  the 
characteristic  at  p.  By  the  smoothness  of  u  there  exists  8^  such  that  each  point  p  +  t.  A,  where 
-hp^f.^bp,  lies  on  a  characteristic.  Let  (a,,  r,,  xj  denote  such  a  characteristic.  Then 
(«!,  fi,  Xj)  is  a  rarefaction  wave  if  the  lines  defined  by  the  characteristics  (uj,  fj,  xj  and  («,,  r,,  xj 
intersect  at  some  f  ^  fp-  This  condition  generalizes  the  ooncqn  of  the  familiar  centered 
rarefaction  wave  for  which  r,  =  fj,  x,  =  x^,  to  non-centered  rarefaction  waves  which  appear  in  the 
generalized  Rionann  problem  solutions  and  also  generalizes  to  a  definition  of  rarefaction  waves  in 
two  spatial  dimensions. 


This  scale  invariance  is  one  of  the  essential  features  of  the  Riemann  problem. 


A  rarefaction  fan  is  an  open  set  in  the  r,  x  plane  all  points  of  which  belong  to  rarefaction 
waves. 

A  shock  is  defined  as  a  imooth  curve  of  points  of  jump  in  the  r,  x  plane  at  each  point  of 
which  the  conditions  (1.4)  and  (1.5)  are  satisfied.  A  shock  curve  (wave)  will  be  denoted 
{ui{t),  uXt),  tj,  Xj)  where  Ui{t)  is  the  limit  value  of  the  solution  u  to  the  left  of  the  shock  wave  at 
time  r  >  r^,  «,  is  the  conesponding  right  limit,  and  (tj,  xj  is  the  origin  point  of  the  shock  wave. 
For  the  Riemann  problem  (1.1)  (2.1),  t,  =  Q,x,  =  Xq. 

A  constant  state  is  a  domain  (connected  open  set)  in  the  r,  x  plane  over  which  the  solution  is 
constant. 

Let  /(«)  i  ChR  "R.  Let  {«i,  1*2}  denote  a  dosed  mterval  where 

{''i.«2}-|[„2,„j     if«i>«2.  ^2.3) 

A  rarefaction  interval  off  shall  denote  an  interval  which  satisfies  either 

/"(«)  >  0      if  Ml  <  M2     for  every  u  d  {u^,  u-^  ,  (2.4a) 

or 

/"(«)  <  0      if  Ml  >  uj     for  every  u  i{u^,u^  .  (2.4b) 

Let  f(Mi,  M2;  u):R  -R  denote  the  convex  envelope  function  for  /(m)  restricted  to  the 
interval  {mi,  M2}  of/  such  that 

JuppCT  convex  envelope  from  Mi  to  M2 ,     if  Mi  >  M2  ,_  „ 

cC"!.  "25 ")  ■  llower  convex  envelope  from  Mi  to  M2 ,     if  Mi  <  M2  ^  '  ' 

Let  C(mi,  M2;  h):  /?  -  /?  denote  the  linear  map  restricted  to  the  domain  interval  {mi,  u-^  such 


that 


i.e.  for  all  M  in  {mi,  u-^, 


dC(Mi,M2;M)       /(m2)-/(«i) 


du  M2  -  Ml         ' 

Denote  the  restriction  of  the  function /(h)  to  the  sub-interval  {mj,  u^  by/(Mi,  M2;  m). 
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Theoron  2.1:  (Oldnilc  [5])  Consider  the  Riemann  problem  (1.1),  (2.1).  Assume 
/(h)  i  C^  :R  -^R  such  that  /  has  a  finite  number  of  inflection  points.  Then  the  unique,  entropy 
obeying  solution  cxmnecting  the  states  «,  and  u^  is  determined  solely  from  £{ui,  u/,  u),  as  described 
below.  Further,  the  solution  is  piecewise  smooth,  with  a  single  irregular  point  corresponding  to 
the  discontinuity  at  jcq  in  the  initial  data. 

The  convex  envelope  function  £(u,,u/,u)  uniquely  divides  the  interval  {«;,  wj  into 
rarefaction  intervals  sqjarated  by  intervals  corresponding  to  the  linear  chord  segments  of  C  Each 
rarefaction  interval  corresponds  to  a  rarefaction  fan  in  the  solution.  The  fans  are  centered  on 
t  =  0,  X  =  xq  and  each  rarefaction  wave  has  a  unique  value  u  corresponding  to  a  state  in  a 
rarefaction  interval.  The  slope  of  the  rarefaction  wave  (v,  0,  xq)  is  given  by  /'(v).  Each  linear 
chord  in  €  covering  an  interval  {u^,  u;^}  corresponds  to  a  shock  wave  (u^,  U2,  0,  xq).  The  solution 
is  bordered  on  the  left  by  the  single  constant  state  u  =  u,  and  on  the  right  by  u  =  u^.  No 
constant  states  appear  within  the  solution  oicept  for  these  two  bounding  constant  states,  which  are 
the  states  u,  and  u^  of  the  initial  data.  The  speeds  of  the  non-linear  waves  found  in  the  solution 
increase  continuously  from  left  to  right  along  the  x  direction.  Fig.  2.1a  shows  an  upper  convex 
envelope  drawn  between  two  states  for  the  function  /(u)  illustrated.  Fig.  2.1b  illustrates  the 
corresponding  Riemann  solution  in  t,x  space,  with  the  three  shocks  (in  this  case)  represented  as 
dark  lines. 

A  solution  given  by  Theorem  2.1  shall  be  referred  to,  for  brevity,  as  the  Oleuuk  solution  to 
{«;,  u^  centered  on  Xq. 

2.1.  The  Generalized  Rkmann  ProMein  for  Initial  Data  With  a  Single  Point  of  Discontiniilty 

Let  A(x,,  xj,  u-^,  U2,  x)  :  R  ->  R  \yt  a  continuous,  differentiable,  1-1  mapping  from  the  interval 
[x^,  X2]  of  the  real  axis  onto  a  rarefaction  interval  {u^,  u-^  of  /  such  that 

A(xi,X2,  Mj,  M2;-«fi)  =  «i  .  07) 

h{x^,X2,U^,U2,X-^  =  "2-  ^  •  ) 

Further,  assume  that  the  derivative  of  li  is  boimded  and  bounded  away  from  0  .  Then  its  inverse 
map  exists,  is  1-1,  onto,  and  has  derivatives  bounded  away  from  0.  Denote  the  inverse  map  by 
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Where  no  ambiguity  exists  the  notation  for  these  maps  shall  be  shortened  to  h(x)  and  A""  («)  . 

Let  /(m)  6  C*:  /?  -  /?  such  that  /"(«)  =  0  at  a  finite  number  of  points.   Let  {a,  b}  and  {c,  d} 
be  rarefaction  intervals  of /(m). 

Definition  2.2:  The  Cauchy  problem 

MliJl  +  3/("(^;  0)  =0  r2  8^ 

dt  dx  '  ^  '  ' 

with  initial  data 


«(x)  = 


a  x-^x,  , 

fti(x„  Xfl,  a,b;x)      X;  <  x  <  ^o  , 
hXxo,  x„  c,  d;x)      Xq<x<x^  , 


(2.9) 


will  be  denoted  as  a  Cauchy  problem  with  generalized  Riemann  data  (generalized  Riemann 
problem,  for  short).  The  unique  solution  to  (2.8)  and  (2.9)  obeying  (1.4)  and  (1.5)  shall  be 
referred  to  as  the  generalized  Oleinik  solution. 

Theorem  2.3:  Let  /(«)  i  C'^.R  ^R  such  that/"(«)  =  0  for  at  most  two  values  of  u.  Then 
the  generalized  Oleinik  solution  to  (2.8),  (2.9)  is  piecewise  smooth;  consists  of  two  constant  states, 
rarefaction  fans,  at  most  three  shock  waves'  and  at  most  two  irregular  points.  A  rarefaction  wave 
either  originates  at  r  =  0  on  the  interval  [x;,  x^]  or  tangentially  to  a  shock  line  at  f  >  0  and  either 
terminates  on  a  shock  line  or  propagates  undistiirbed  for  all  r.  At  most  two  points  can  be 
irregular,  the  initial  discontinuity  in  the  data  at  (f  =  0,  x  =  Xq)  or  the  meeting  pwint  of  three 
shocks'. 

2.1.1.  Proof  of  Theorem  2.3 

The  following  definitions  will  prove  useful.  Let  {mj,  M2}  be  a  rarefaction  interval  of  /  and 
/»(x],  X2,  «i,  U2,x)  a  map  as  in  (2.7).    At  each  point  x  €  [xj,  xj  construct  the  rarefaction  wave 

(u,  0,  x)  with 

'  By  definition  of  a  shock  wave  as  a  smooth  curve,  the  occurences  of  either  two  shock  lines  interacting  to 
produce  a  single  shock  curve,  or  a  single  shock  curve  splitting  into  two  shock  curves,  are  each  described  as 
involving  three  shock  waves. 
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x(r)  =  or  + /'(/.(x))  r  . 

The  set  of  waves  thus  constructed  on  [x^,  xj  form  a  rarefaction  fan  which  shall  be  referred  to  as 
the  rarefaction  fan  determined  by  h{x^,  x-^,  u^,  u^,  x). 

Let  {mj,  «2}  be  a  rarefaction  interval  of/.  For  every  u  in  {mj,  «2}  construct  the  rarefaction 
wave  (m,  0,  Xq)  with 

The  set  of  waves  thus  constructed  form  a  rarefaction  fan  which  shall  be  referred  to  as  determined 
at  Xq  by  the  interval  {u^,  u^. 

Let  {a,  fc}  be  a  rarefaction  interval  off.  Let  A(po,  pj,  a,  £>;  p)  be  a  map  as  in  (2.7)  from  the 
real  interval  [po,  pj]  to  {a,  b}.  Qsnstruct  a  curve  segment  (r(p),  x(p))  from  (t^,  Xq)  to  (t^,  x{)  in 
the  r,  X  plane  such  that 

(^0.  -to)  ■  ('(Po).  Jc(po))  , 

(rj,  xi)  -  (r(pi),  x(pi)) , 

dr(pVdr(p)=/'(/.(po,Pi,a,*;p)). 

Al  each  point  (f(p),  x(p))  construct  the  rarefaction  wave  (/i(p),  t(p),  x(p)).  Then  the  set  of  waves 
constructed  on  the  curve  segment  form  a  rarefaction  fan  which  shall  be  referred  to  as  the 
rarefaction  fan  determined  by  the  curve  (t(p),  x(p)j. 

Lemma  2.4:  Let  >/  be  a  rarefaction  fan  in  the  r,  x  plane  for  some  solution  u(t,  x).  Let 
[x^,  xj  be  an  interval  in  M  parallel  to  the  x-axis  at  some  time  t  such  that  the  rarefaction  wave 
passing  through  the  point  (f,  x^)  has  value  u  =  u^  and  the  rarefaction  wave  passing  through  the 
point  (f,  X2)  has  value  u  =  M2-  Then  {mj,  u^  a  a  rarefaction  interval  of  /  and  the  fan  defines  a 
continuous,  1-1  mapping  h'  (xj,  X2,  Uj,  U2;x)  -.R  -R  from  [xj,  xj  onto  {u,,  u-^  whose  derivatives 
are  bounded  and  bounded  away  from  0  . 

The  proof  for  Lemma  2.4  follows  from  the  definition  of  a  rarefaction  wave.  □ 

The  proof  for  Theorem  2.3  is  constructive.  All  possible  Cauchy  data  (2.9)  fall  into  four 
classes: 
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Cl:    {a,  b}  and  {c,  d}  are  not  disjoint,    /"  >  0  on  {a,  b}  U  {c,  </}  , 

C2:    {a,  ft}  and  {c,  </}  are  not  disjoint,    f"  <0 on{a,  b}  U  {c,  d}  , 

C3:    {a,  ft}  and  {c,  i}  are  disjoint,    ft  >  c  , 

C4:    {a,  ft}  and  {c,  d}  are  disjoint,    ft  <  c  . 

Qasses  C3  and  C4  can  each  be  divided  into  four  subgroups: 


51 

a<b. 

c  <d  , 

52 

a>b, 

c>d  , 

53 

a  <b. 

c>d, 

54 

a>b. 

c<d. 

For  classes  CI,  C2  there  are  no  inflection  points  of  /  in  the  interval  [a,d].  For  subgroups  SI  and 

52  of  classes  C3  and  C4  there  are  either  0  or  2  inflection  points  in  the  interval  [a,d].  Subgroups 

53  and  S4  of  classes  C3  and  C4  contain  one  inflection  point  in  the  interval  [a,  d].  Theorem  2.3 
will  be  proven  for  the  case  of  0  and  1  inflection  point  in  the  interval  [a,  d]  of/.  The  proof  for  the 
case  of  two  inflection  points  will  be  omitted  for  brevity. 

/  has  Zero  Inflection  Pofaits 

Restriction  of  /  to  no  inflection  points  implies  initial  data  from  either  classes  CI,  C2  or 
subgroi^K  SI  or  S2  of  C3  and  C4.  We  study  data  of  classes  CI,  C3.S1  and  C3.S2;  the  proofs  for 
classes  C2,  C4.S1  and  C4.S2  follow  analogously. 

QassCl 

There  are  four  arrangements  of  initial  data  in  the  dass  CI: 

Cl.Al  a  <  c  <  b  <  d  , 

C1.A2  a  <  c  <  d  <  b  , 

C1.A3  c  <  a  <  d  <  b  , 

C1.A4  c  <  a  <  b  <  d  . 

The  generalized  Oeinik  solution  to  Cauchy  data  of  type  Cl.Al  is  examined  below. 

Let  {a,  ft}  and  {c,  d}  be  rarefaction  intervals  of  type  Cl.Al  .  Consider  maps 
hi(xi,  Xq,  a,  ft;  at)  and  H^xq,  x^,  c,  d;x)  as  in  (2.7).  Construct  the  rarefaction  fans  determined  by 
hi(x)  and  h^x).  Then  through  each  point  p  -  (r,  x)  in  the  plane  such  that 

xo  +  /'(c)  t:SxrSxo+  fib)  t 
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there  passes  a  unique  rarefaction  wave  originating  from  the  interval  [x;,  xj  and  a  unique 
rarefaction  wave  orig  nating  from  the  interval  [xq,  xJ,  thus  defining  respectively  the  unique  values 
Ui(p)  and  u^(p). 

Consider  the  curve  (f,  x(f))  defined  implicitly  by 

^('^  =  ^0  +  /J  ^  'i'  ,  (2.10) 

where 

dx  _Mit))-fiuM  ._,. 

dt         uM-uXt)     '  ^^"^ 

with  Ui{t)  and  u^t)  determined  from 

x(r)  =  h,-^  ix„  xo,  a,  b;  u^t))  +  fiu^it))  t ,  (2.12) 

x(0  =  \-i  (xo,  x„  c,  d;  uXt))  +  /'(«.(0)  '  •  (2.13) 

We  remark  that  for  {a,  b)  and  {c,  d\  of  type  Cl.Al,  the  set  of  equations  (2.10)  through  (2.13) 
are  well  defined  for  all  t  and  the  curve  (r,  x(r))  is  smooth. 

Proposition  2.5:    For  the  curve  {t,  x(f))  defined  by  (2.10)  through  (2.13): 

d  u,(t)  d  uXt) 

1)  '^'  <  0  ,   — f^  >  0,  Vf  >  0  . 
dJ  dt 

2)  There  exists  g  €  [c,  b]  such  that  U;(f)  -  g  ,  m^O  -  S   as  f  -  <». 

3)  The  curve  defined  in  (2.10)  through  (2.13)  considered  as  a  curve  of  points  of  jump  u,(t)  -  uXt) 
obeys  (1.4)  and  (1.5). 

Proof:    Note 

/'    u,(f))  >  0,      (hrKu))'  >  0,     V«,(r)  on  [a,  b]  , 
f"iuM>0,     (V'(«))'>0,     yfuXi)  on  [c,d]. 
Equations  (2.12)  and  (2.13)  yield 

If'iuXt))  -  f(uM]  t  =  K-\uXt))  -  hrKuM)  .  (2.14) 

For  r  >  0,  the  positivity  of  the  right  hand  side  of  (2.14)  implies 

/'(«/W)-/'(«.(0)>0.  (2.15) 
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By  the  requirement  that  {a,  b}  and  {c,  d}he  of  class  Cl.Al,  (2.15)  implies  (using  2.11) 

/'(«/(0)>-f  >/'(«r(0)-  (2.16) 

Therefore,  the  time  derivatives  of  (2.12)  and  (2.13), 

^  =  KhrKu,))'  +  /"(«;)  t]  «/(f)  +  /'(«,)  ,  (2.17) 

^  =  [('«r^(«.))'  +  /"(«.)  f]  «/(0  +  /'(".)  ,  (2.18) 

lead  to  the  conclusion 

«/(f)  <  0  ,  «/(0  >  0  ,  Vr .  (2.19) 

Further  (2.16)  and  (2.19)  imply  the  existence  of  g  €  [c,  b]  such  that  as  f  -« oo, 

uM-g,    uXt)-g,    ^-f'is). 

By  construction,  the  curve  {t,  x{t))  considered  as  a  curve  of  points  of  jump  obeys  (1.4).  It  is 
easy  to  verify  that  the  entropy  condition  (1.5)  is  also  satisfied.  From  Proposition  2.5  note  that  the 
strength  of  the  shock  decreases  monotonically  to  zero  as  r  -  ao.  n 

All  rarefaction  waves  originating  in  the  interval 

for  the  construction  in  Proposition  2.5  intersect  the  curve  (f,  x(t)).   Terminate  these  waves  at  the 
curve.   All  rarefaction  waves  originating  in  the  intervals 

A;- 1(a)  5  X  ^  hrKg) 
h-Kg)  ^x^  h-\d) 

propagate  unimpeded  for  all  r  >  0.    The  rarefaction  fans  and  shock  curve  thus  defined  by  this 
construction  shall  be  said  to  be  determined  by  the  functions  h^x,,  Xq,  a,  b;  x)  and  h^x^,  x^,  c,  d;  x). 

Proposition  2.6:  Let  {a,  b}  and  {c,  J}  be  as  in  Cl.Al  .  Then  the  generalized  Oleinik 
solution  u(f ,  x)  to  (2.8)  (2.9)  is  given  by 
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the  constant  state  u  =  a  for  x  ^  X)  +  f'(a)  t  ,  r  &  0  , 

the     shock     curve     and     rarefaction     fans     determined     by     A;(x)     and     h^x)     for 
X,  +  f'ia)  t<x<x,  +  fid)  t  ,     r  a:  0  , 

the  constant  state  u  =  rf  for  x^  +  f'(d)  t  ^  x  ,  f  &  0  , 

The  proof  follows  immediately  from  Proposition  2.5.  □ 
A  sketch  of  this  solution  is  given  in  Fig.  2.2  . 

The  generalized  Oleinik  solution  to  Cauchy  data  of  types  C1.A2  through  C1.A4  have  exactly 
the  same  structure  as  for  Cl.Al  .  For  each,  the  strength  of  the  single  shock  curve  (t,  x{t))  in  each 
solution  decreases  raonotonically  to  zero  as  f  -  <»,  and  dx  /  dt  "  f'{g)  where  for 

C1.A2  c  <  g  :s  d  , 
C1.A3  a  <  g  s  d  , 
C1.A4      a  <  g  :s  b  . 

Qass  C3.S1  -  No  Inflection  Point 

Let  {a,  b)  and  [c,  d}  be  rarefaction  intervals  of  type  C3.S1  such  that  for  every 
V  €  {a,  b},  w  ^{c,d} 

f(v,  w;u)  =  C(v,  w;u)  . 

Then  /'(v)  >  f'(w)  for  every  choice  v,  w.  Consider  the  curve  (f,  x(f))  inyjlidtly  defined  by 
(2.10)  through  (2.13).  Statement  1  of  Proposition  2.5  follows  from  (2.14)  through  (2.19). 
However,  the  disjointness  of  the  intervals  implies  there  exists  finite  times  //  >  0,  t^>  0  such  that 
u,(t,)  -  a  and  M^(fJ  -  d.  For  all  f  >  f;  let  «,(')  "  '^  ^^  for  all  r  >  f^  let  u^t)  «  d.  With  this 
extension,  (2.10)  through  (2.13)  are  defined  for  all  t.  The  curve  (r,  x{t))  thus  defined,  considered 
as  a  curve  of  points  of  jui  p  u,(t)  -  «^(f)  obeys  the  conditions  (1.4)  and  (1.5).  In  contrast  to  class 
CI  data,  the  shock  strength  attains  a  constant  value,  |a  -  d|,  as  r  -  ». 

All  rarefaction  waves  originating  from  the  interval  [x/,  x^]  intersect  the  shock  curve  and  are 
terminated  there.  The  shock  curve  and  rarefaction  fans  so  defined  (together  with  the  bounding 
constant  states  u  —  a,  u  =  d  will  also  be  referred  to  as  determined  by  the  functions  h^{x)  and  h^(x) 
since  their  definition  is  merely  an  extension  of  the  usage  defined  after  Proposition  2.5. 
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Proposition  2.7:  Let  {a,  b}  and  {c,  J}  be  disjoint  rarefaction  intervals  on/  such  that  for  all 
V  £  {a,  b}  ,  w  (.  {c,  d} 

e(y,w;u)  =  C(y,w;u)  . 

Then  the  generalized  Oleinik  solution  to  (2.8)  (2.9)  is  given  by  constructing  the  shock  curve  and 
rarefaction  fans  determined  by  /»j(x;,  Xq,  a,  b;x)  and  hX^Q,  x^,  c,  d;  x). 

Proof:  The  above  construction  proves  this  proposition  for  rarefaction  intervals  in  class 
C3.S1,  vdiere  the  condition  /"  >  0  was  used  implicitly  in  (2.17)  and  (2.18).  With  trivial 
modifications,  analogous  arguments  can  be  used  to  prove  this  proposition  for  appropriate  data  in 
aasses  C3  or  C4. 

Qass  C3.S2  -  No  Inflection  Point 

Data  of  class  C3.S2  insures /'(fc)  s  /'(c). 

Proposition  2.8:  Let  {a,  b]  and  {c,  </}  be  rarefaction  intervals  of  /  of  class  C3.  Further 
assume /'(fc)  s  /'(c).  Then  the  generalized  Oleinik  solution  u{t,  x)  to  (2.8)  (2.9)  is 

The  constant  state  u{t,  x)  -  a  iox  x  <  Xf  +  f'(a)  t,  f  >  0  , 

The  rarefaction  wave  determined  by  A^(x;,  xq,  a,  b;x)  for 

X,  +  f'{a)  r  <  X  <  xo  +  fib)  t  ,  r  &  0  , 

The  Oleinik  solution  to  {b,  c}  centered  on  Xq  for  Xq  +  /'(*)  t  <  x  <  Xq  +  f'(c)  t  ,  r  &  0  , 

The  rarefaction  wave  determined  by  h^xo,  x^,  c,  d;x)  for 

xo  +  /'(c)  r  <  X  <  X,  -t-  f'(d)  t  ,  t^O, 

The  constant  state  u(f,  x)  =  <f  for  x  <  x^  +  f'(d)  t,  f  &  0  . 

Proof:    The  proof  is  trivial;    the  condition  f'(b)  ^/'(c)  guarantees  that  the  rarefaction  waves 
determined  by  A;  and  h^  do  not  interact. 

Thus  all  solutions  for  allowed  Caudiy  data  with  /  restricted  to  have  no  inflection  point  are 
pnecewise   smooth,   with   at  most   a  single  shock  and  at  most   a  single  irregular  point   at 

(f    =    0,  X    =   Xq). 
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/  has  One  Inflection  Point 

Restriction  of  /  to  a  single  inflection  point  implies  initial  data  from  subgroups  S3  or  S4  of  C3 
and  C4.  We  study  data  of  C3.S3  and  C3.S4;  the  proofs  for  C4.S3  and  C4.S4  follow  analogously. 

aassC3.S3 

Lemma  2.9:    Let  /,  {a,  b},  and  {b',  a'}  be  chosen  such  that:  {a,  fc}  is  a  rarefaction  interval 
of/ on  which /"(«)>  0,  {i',  o'}  is  a  rarefaction  interval  of/ on  which /"(«)<  0,  a'<b,and 

e{b,  b';  u)  =  C(*,  b';u),  fib')  =  ^<^^^'*'i")  , 

du 

£(a,  a';u)  =  C(a,  a';  u)  ,  f'(a')  =  ^C^a^''^")  . 

For  every  v  €  {a,  b}  there  exists  vi'(v)  6  {b',  a'}  such  that 

av,  »v(v);  u)  =  Civ,  »v(v);  u)  ,  /'(>v(v))  =  ^  ^(^  *^(^);  ")  . 

du 

Then 

1)      K'(v)  is  a  1-1  map  from  {a,  b}  onto  {b',  a'}, 

2)  /-(^,))j^^/'Kv))-r(v) 

dv  ^(V)  ~  V 

Proof: 

1)  This  follows  from  the  definition  of  h'(v). 

2)  For  all  V  €  (a,  A),  8v  sufficiently  small, 

/(»v(v))  -  /(v)  =  /'(h-(v))  [,v(v)  -  v]  ,  (2.20) 

f{yv{\  +  8v))  -  /(v  +  8v)  =  /'(^(v  +  8v))  [^(v  +  8v)  -  v  -  6v]  ,  (2.21) 

Subtracting  (2.20)  from  (2.21),  arranging  terms  and  taking  8v  -  0 

/.(,(,))   ^    _  ^'(v)    =  /"(W(V))   ^   »V(V)   -H  /'(H.(V))   ^    -  /"(H.(V))   ^   V   -  /'(»V(V))   , 

from  which  the  conclusion  2)  follows.  D 

Let  {a,  b}  and  {fe*,  a*}  be  as  in  Lemma  2.9.    Let  h^{xi,  xg,  a,  i;x)  be  the  usual  map  from 
the  rarefaction  interval  {a,  b\  to  \x^,  xj.   Consider  the  curve  (f,  Jt(f))  defined  by 


16 


dx 


'(0  =  ^0  +  Jo  -f  *  .  (210) 

^=/'(h'(«XO)).  (2.22) 

with  H'(uj(f))  defined  in  Lemma  2.9  and  U;(0  defined  by 

x{t)  =  Ar'  (^i.  ^0.  «.  *;  «X0)  +  /'(«;(0)  ^  •  (2.12) 

There  exists  a  finite  time  f y  >  0  at  ^^cfa  u^{t^  =  a  and 

.|=/'Ka))=/V). 

For  every  t  >  t,  let  u^t)  =  a.  Then  (2.10),  (2.12)  and  (2.22)  are  defined  for  all  t. 

Proposition  2.10:    The  curve  defined  by  (2.10),  (2.12)  and  (2.22)  cxmsidered  as  a  curve  of 
points  of  jump  Uj  -  w(u,)  obeys  (1.4)  and  (1.5). 
Proof:    Follows  from  the  construction,  n 

Construct  the  rarefaction  fan  determined  by  A;.  If  a  wave  from  this  fan  intersects  the  above 
shock  curve,  terminate  the  wave  on  the  curve.  For  a  point  (r^,  x(f^))  on  this  shock  curve  where 
0  ^  f^  s  t„  there  exists  a  uniquely  defined  uXQ  =  ^("lih)-  ■^so  note  from  (2.22)  and  Lemma 
2.9 

x"it)  =  f'Xu,))  »v'(«/)  «/(0  =  ^'^''jrfllf^"'^  «/'(')  <  0-  (2-23) 

The  rarefaction  fan  determined  by  this  curve  can  therefore  be  constructed.   Note  further  that  the 
'left  bounding'  wave  (a',  t,  x{t))  of  the  fan  propagates  parallel  to  the  shock  curve. 

The  two  rarefaction  fans  and  shock  wave,  together  with  the  bounding  constant  state  u  =  a, 
thus  constructed  from  (2.10),  (2.12),  (2.22)  and  (2.23)  shall  be  denoted  as  determined  by  the 
Junction  hjix)  and  the  state  u  =  b' . 

Let  /  have  a  single  inflection  point  and  {a,  b}  and  {c,  J}  be  rarefaction  intervals  of  type 
C3.S3  .  Then  there  exists  a  rarefaction  interval  {b',  a'}  as  described  in  Lemma  2.9.  Three 
possibilities  exist  for  the  orientation  of  {c,  </}  with  respect  to  {b',  a'}  : 

Case  1:    {c,  d\  and  {b',  a'}  are  disjoint  with  d  ^  a'  . 

Case  2:    {c,  d}  and  {b',  a'}  are  disjoint  with  c  ^  b'  . 
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Case  3:    {c,  </}  and  {fc*,  a'}  are  not  disjoint 
For  the  case  J  &  a*,  the  generalized  OleLiik  solution  u(t,  x)  to  (2.8)  (2.9)  is  given  by  Proposition 
2.7. 

For  the  case  c  -^  b' ,  the  generalized  Oleinik  solution  u(t,x)  to  (2.8)  (2.9)  is  given  by  the 
following: 

The  constant  state  u  =  a,  two  rarefaction  fans  and  the  shock  curve  determined  by  the 
function  h^x)  and  the  state  u  =  b'  for  x  ^  f'{b')  t  ,  r  >  0  . 

The  rarefaction  wave  determined  at  x^  by  the  interval  {c,  b')  for  /'(£»*)  t  ■^  x  ^  f'(c)  t  , 
t>0. 

The  rarefaction  wave  determined  by  the  h^x)  for  f'(c)  t  ^  x  ^  f'(d)  t  ,  t  >  0  . 

The  constant  state  u  =  d  for /'(c)  t  ^  x  ,  t  >  0  . 
This  solution  is  sketched  in  Fig.  2.3  . 

For  the  case  in  which  {c,  d}  and  {b',  a'}  are  not  disjoint  we  remark  that  the  generalized 
Oleinik  solution  can  be  described  by  a  solution  similar  to  case  1  above  f or  0  ^  f  ^  r^  and  by  a 
solution  similar  to  case  2  for  f  <  f.  The  time  t  is  determined  by  the  propagation  of  a  rarefaction 
wave  for  some  state  u  -  g  where  g  €  {i",  a'}  D  {c,  d}. 

aassC3.S4 

Lemma  2.9  and  Proposition  2.10  can  be  written  analogously  for  data  of  type  C3.S4,  with  the 
arguments  revolving  around  an  interval  {c',  cT].  Again  the  construction  of  a  solution  to  the 
generalized  Riemann  problem  can  be  broken  into  three  cases  and  dealt  with  as  above. 

The  solutions  for  classes  C4.S3  and  C4.S4  follow  in  the  same  manner.  Thus  we  have  shown 
that  the  solutions  to  the  generalized  Riemann  problem  for  /  having  at  most  a  single  inflection 
point  are  piecewise  smooth  and  consist  of  at  most  a  single  shock  wave  and  at  most  a  single 
irregular  point  (at  jq). 

/  has  Two  Inflection  Points 

In  the  interest  of  brevity,  we  omit  the  solution  for  the  case  of  /  having  two  inflection  points. 
The  method  of  proof  is  repetitive  to  that  presented  above.   The  proof  shows  that  in  addition  to 


-18- 

solutions  with  a  single  shock,  either  two  shock  waves  can  be  present  for  all  r  >  0  or  three  shock 
waves  can  appear  in  the  solution  [i.e.  a  single  shock  which  is  present  f or  0  ^  r  ^  r^  which  then 
'splits'  (via  interaction  with  rarefaction  waves)  into  two  shocks  \^cfa  are  present  for  t,  ^  t]. 

2.2.  A  Coi^ecture  on  Loss  of  Piecewiae  Smoothness  for  More  Than  Two  Inflection  Points 

The  constructive  proof  given  for  Theorem  2.3  is  not  readily  extendable  to  /  having  three 
inflection  points  since  the  number  of  cases  to  be  considered  is  large.  In  fact,  we  argue  that  for 
/(m)  having  three  or  more  inflection  points,  the  solution  need  not  be~piecewise  smooth.  An 
example  is  sketched  below  in  wiiicfa  a  solution  having  a  countably  infinite  number  of  arbitrarily 
small  smooth  pieces  may  occur. 

Let  {a,  b}  and  {c,  </}  be  rarefaction  intervals  of /(«)  as  illustrated  in  the  lower  left  hand  side 
of  Rg.  2.4  .  /(u)  has  three  inflection  points  in  the  interval  [d,  b].  Consider  the  corresponding 
generalized  Riemann  problem  (2.8),  (2.9).  Using  the  construction  methods  of  the  proof  of 
Theorem  2.3,  it  is  seen  that  the  generalized  Oleinik  solution  initially  has  a  single  shock  curve 
determined  by  hi(x)  and  hXx)  starting  at  Xq  with  slope  d  C(b,  c;  u)  I  du.  The  rates  at  which  the 
left  and  right  states  «;(r)  and  uXt)  change  along  the  shock  curve  depend  on  the  functions  hi{x)  and 
h^x).  Given  {c,  d\,  (and  hence  {a,  b}),  [xq,  ocJ  and  h^x),  it  is  possible  to  determine  A^x)  and 
[x;,  xj  such  that  for  all  time  t  >  0 

eiuXt),  uXt);  u)  =  C(u,it),  git);  u)  U  edit),  uXt);  u)  ,  (2.24) 

where 

,u  rA^       dCiuXt),8it);u)  _  dCigit),uXt);u) 

f  iiit)) Jj J^ (2.25) 

The  convex  envelope  (2.24),  though  consisting  of  two  diords,  defines  a  single  shock  line  and 
represents  a  metastable  balance  between  a  single  and  a  two  shock  configuration.  This  construction 
can  be  perturbed  to  aeate  a  single  'bubble';  that  is  a  time  interval  fj  <  f  <  f,  over  which  the 
shock  separates  into  two  shocks  as  is  illustrated  on  the  right  in  Fig.  2.4  .  This  is  accomplished  by 
adjustment  of  the  rate  at  which  the  rarefaction  waves  (i.e.  adjustment  of  the  respective  h 
functions)  enter  the  shock  waves  from  the  extreme  left  and  right.  The  above  construction  can  be 
iterated  to   produce   a  second  bubble  starting  at  t-^  by  adjusting  (slowing)   the  rate  A/(x) 
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appropriately.  It  has  not  as  yet  been  determined  if,  for  a  given  set  of  intervals  {a,  b},  {c,  d}, 
[xj,  xq],  and  [xq,  jJ,  is  possible  to  construct  a  sequence  of  generalized  Riemann  problems,  where 
the  j'th  problem  is  characterized  by  the  choice  of  functions  h'i(x),  h'Xx),  for  which  the  solution  to 
the  j'th  problem  has  i  such  bubbles.  If  this  was  possible,  the  solution  to  the  limiting  problem 
obtained  as  /  -  oo  would  not  be  piecewise  smooth.  We  note  it  is  certainly  possible  to  create  an 
infinite  sequence  of  such  problems  if  only  {a,b},  {c,  </},  and  [x^,  xj  are  held  constant  and  the 
interval  [xq,  xJ  is  allowed  to  vary  with  i.  It  remains  to  be  proven  that  such  a  sequence  of 
problems  can  be  diosen  such  that  [xq,  xJ  also  remains  fixed. 

3.  The  Generalized  Riemann  Problem  for  Initial  Data  with  Multiple  Points  of  DiscontiDiiity 

Let  X;  <  Xq  <  Xi  <   ■  •  •   <  x„  <  x^ .  Consider  (2.8)  with  initial  data 


"(0,  x)  = 


U,  X  <Xi  , 

h,ix„  X(j,   «;,   Uq,  x)  X;    <  X    <   Xq   , 

Ay(Xy,Xy+i,   «y,   «y+i;x)  Xy    <  X    <  Xy+ j    ,     >    =    0,    1 ,      ••      ,«"!,  (3.1) 

hXx„,  x„  u„,  u/,  x)        x„  <  X  <  X,  , 
u_  x,^  X  , 


where  hj{x),  hj(x),  hXx)  are  as  in  (2.7),  and  {u„  Uq},  {uq,  uJ,  ....  {u„,  u^}  are  rarefaction  intervals 
of  /.  The  Cauchy  problem  (2.8),  (3.1)  is  an  extension  of  the  definition  of  the  generalized 
Riemann  problem  to  include  multiple  initial  discontinuities. 

Theorem  3.1:  The  solution  to  (2.8)  with  initial  data  (3.1)  (or  f  d  C^  -.R  ^  R  having  at  most 
one  inflection  point  is  piecewise  smooth.  There  are  no  more  than  (n  +  l)I  shocks  in  the  solution 
and  no  more  than  (n  +  1)!  irregular  points. 

The  restriction  to  either  zero  or  one  inflection  point  for  f(u)  guarantees  that  when  two 
shocks  interact  a  single  shock  is  produced.  Thus  given  a  finite  number  m  of  initial  discontinuities, 
the  number  of  shocks  ii  Jie  solution  remains  bounded  by  m!  and  the  solution  will  be  piecewise 
smooth.  For  /  having  two  or  more  inflection  points  there  is  no  a  priori  bound  on  the  growth  of 
the  number  of  shocks  in  time  since  an  interaction  of  two  (incoming)  shocks  can  conceivably  give 
rise  to  two  or  more  (outgoing)  shocks  (or  to  a  single  shock  which  later  splits  via  interaction  with 
rarefaction  waves  into  two  or  more  shocks)  and  the  possibility  of  generating  a  solution  that  is 
richer  in  structure  than  piecewise  smooth  exists.  It  is  only  through  the  decease  in  variation  of  the 
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solution  (ie.  the  decrease  in  shock  strength  with  in  time)  that  such  a  bound  may  be  possible.  Such 
a  bound  on  the  growth  rate  of  shock  number  combined  with  a  general  statement  for  the  >  2 
inflection  point  case  of  §2  is  required  for  a  more  general  statement  on  the  piecewise  smoothness 
of  the  solution  to  the  generalized  Riemaim  problran  for  initial  data  with  multiple  points  of 
discontinuity. 

Theorem  3.1  will  be  proven  by  finite  induction  on  n  in  (3.1).  The  case  n  =  0  is  the 
generalized  Riemann  of  §2.1  .  The  case  n  =  1  is  given  below.  It  suffices  to  illustrate  the  salient 
points  of  the  general  induction  step. 

Proof  of  Theorem  3.1  for  n  =  1:  Consider  equation  (2.8)  under  three  possible  cases  of 
Cauchy  data: 


u(x)^ 


h,(xi,Xo,u,Uo,x) 


X  ^x,  , 

Xi<X  <Xo, 


<x  <x 

x,sx, 


1  » 


(3.2) 


u(x)  = 


«0 


ho(.Xo,x^,U(t  u^;x) 
A,(xi,x„«i  «,;x) 


Xq  <  a:  <  or 

X,    <  X  <  X, 


1  > 


(3.3) 


«(x)  = 


hi(xi,Xo,UiUo,x) 
/»o(xo,xi,  Uo_«i;x) 
hXxi,x„  u^u/.x) 


X;    <  X   <   X, 

Xfl  <  X  <  X 

Xi  <  X  <  x^ 
X,  S  X  . 


0  » 

1  > 


(3.4) 


From  §2,  construct  the  solutions  to  (3.2)  and  (3.3),  denoting  them  ui^t,  x)  and  u^it,  x) 
respectively.  Denote  the  solution  to  (3.4)  as  «c(f,  x).  The  following  four  (exhaustive  but  not 
mutually  exclusive)  possibilities  exist  for  the  solutions  u^  and  Uji.* 


*  As  xq  and  xi  are  spatially  separated,  intuitively  uc  is  expeaed  to  be  campcsed  thusly: 

uc  =  left  piece  of  «/,  U  piece  due  interaction  of  waves  xq  and  xj  U  right  piece  of  u^  . 
Conditic3os  PI  through  P4  characterize  the  interactions  between  waves  from  xq  and  xj  . 
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PI)    There  exists  a  rarefaction  wave  (v,  0,  x^,)  common  to  both  solutions  with  v  €  {uq,  mJ, 
x^  i  (xq,  Xj)  which  propagates  unimpeded  for  all  t. 

P2)    Both  M^  and  u^  have  a  shock  curve  and  there  exists  a  point  pj  =  ('2'  ^2)  common  to  the 
shock  curve  of  each  solution  for  some  f2  >  0  . 

P3)    u^  has  a  shock  curve,  denoted  (u/(0,  u}(t),  x)  and  m^C's)  =  «i    for  some  f3  >  0  . 
P4)    u^  has  a  shock  curve,  denoted  {uf{t),  u;(f),  x)  and  u^it^)  =  Mq    for  some  fj  >  0  . 

If  PI  holds  (if  neither  solution  has  a  shock  wave,  then  condition  PI  holds)  then 


.    r,  ,^  _  /  "i(^'  t)  ii   x^x,+  /'(v)  f  , 
"cC^-'^      [  ««(x,0  if   x>x,+/'(v)r. 


(3.5) 


is  the  unique  solution  to  (2.8)  (3.1)  obeying  (1.4)  and  (1.5). 

If  PI  does  not  hold  let  f^^^  denote  the  minimum  of  the  existing  values  r,,  fj,  14  . 

K  f^„  =  u  then  let  ^l  -  (ufit),  «0(r),  0,  Xq)  and  2^  -  (u}(t),  uHt),  0,  x^)  denote  the  shock 
waves  propagating  respectively  from  xq  in  u^  and  Xj  in  m/j.  Then  for  f  s  fj 


"^("'  '^      1  «,(x,  r)  if 


X^V'(V2)+/'(V2)'. 


(3.6) 


X>  V'(v2)+/'(v2)r, 

where  Vj  =  «r('2)  ~  "/('i)  '^  "^  ^^  interval  {uq,  u^}.    For  r  >  r,.  wc(^>  0  '*  *^  generalized 
Oleinik  solution  to  (2.8)  with  initial  data  (at  12)  __ 


«('2.  Jf)  = 


X  S  X 


/    ' 


where 


a/'  (Xj,  x/S  ^^^(fj),  U/,  X)       Xj  <  X  <  x/'  , 


x/'  =  X,  +  /'(u,)  fj  . 


A/Xx)  is  defined  by  the  rarefaction  interval  {«;,  u]  (rj}  where  u]{t-^  is  the  left  limit  of  u  at  (t,,  x;) 
for  the  shock  curve  Zq  in  u^  propagating  from  x^  in  (3.6).  (Zq  is  identical  with  X^  for 
0  ^  »  ^  f,  •)  Similarly  h^\x)  is  defined  by  the  rarefaction  interval  {u^,  up  (/J}  where  u^^tj)  is  the 
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right  limit  of  u  at  (fj,  x^  for  the  shock  curve  S^  propagating  from  x^.   (Si  is  identical  with  S^  for 
0  s  f  s  fj  .) 

If  r^j,  =  fj,  let  X3  =  xi  +  /'(mi)  t  .  There  exist  two  possibilities  for  the  solution  uciiu^ha& 
a  shock  wave.  (If  m^  has  no  shock  the  solution  for  u^  is  deduced  »-asily  from  the  discussion 
below.)  These  two  possibilities  are  dther  that  m^  has  a  leftmost  rarefacti   >  fan 

a)  centered  at  the  point  x^  (Fig.  3.1a), 

b)  determined  by  the  shock  wave  S^  originating  at  (0,  x{)  (Fig.  3.1b). 

We  note  that  in  either  case  this  leftmost  rarefaction  fan  in  u^;  is  a  smooth  continuation  of  the 
rightmost  rarefaction  fan  in  the  solution  u^  f or  f  ^  ty  For  t  s  fj,  let 

At  (ty,  X3)  the  shock  Zq  originating  from  (0,  xq)  encounters  the  first  rarefaction  wave  from  the 
solution  for  x  >  xj.  For  case  a)  above,  the  shock  1q  is  determined  for  f  >  fj  by  the  rarefaction 
fan  determined  by  hi(x„  Xq,  u,,  Uq;x)  and  the  rarefaction  fan  centered  on  (0,  x^).  Etenote  the 
curve  followed  by  the  shock  Zq  for  f  >  f,  as  (f,  Xjjj)).  Then,  for  case  a), 


«c('.  ^)  = 


"R    Jc>xj^(r),  3 


In  case  b)  the  shock  2o  is  determined  for  r  >  fj  by  the  rarefaction  fan  determined  by 
hiixj,  Xq,  Ui,  Uq;x)  and  the  rarefaction  fan  determined  by  the  shock  line  Xj.  In  either  case  a)  or 
b),  the  resulting  shock,  1q  hits  the  shock  Xj  at  some  time  fj  >  f 3  or  it  does  not.  If  it  does  not,  the 
solution  f or  all  f  >  f,  is  given  by  (3.8).  If  it  does,  the  solution  f or  fj  s  r  :s  r,  is  given  by  (3.8) 
and  f or  f  &  ^5  is  given  by  the  solution  to  a  generalized  Riemann  problem  of  the  type  discussed  in 
52.1.  It  is  to  be  noted  that  the  solution  is  C^  continuous  at  fj  (ie  all  waves  propagate 
continuously  from  '3  -  t  to  ^3  +  t  and  all  waves  except  the  interacting  shock  lines  are  C^ 
continuous  at  ty 

If  f      =  r4,  an  analogous  construction  to  that  for  case  P3  holds. 
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In  all  cases  the  constructed  solution  u^  is  thus  piecewise  smooth.  The  general  induction  step 
for  Theorem  3.1  follows  the  above  argument  with  a  more  complicated  notation,  n 

4.  The  2-Dbiieiisioiial  Riemann  Problem 

For  convenience  we  shall  use  the  notation /^  —  /,  /j  -  g,  Xj  -  x,  ^2  -  y  for  the  discussion 
in  two  spatial  dimensions. 

Definition  4.1:  The  Caudiy  problem  (1.1)  with  initial  data  piecewise  constant  on  a  finite 
number  of  wedges  focused  on  a  single  point  in  the  x,  y  plane  shall  be  referred  to  as  the  two 
dimensional  Riemann  problem.  Wthout  loss  of  generality,  this  point  can  be  taken  to  be  the  point 
X  =  0,  3^  =  0. 

For  the  case  g  «■  /,  under  the  45°  rotation  2^  =  x  +  y,  2ti  =  y  -  x,  (1.1)  becomes 

4^  +  -^  =  0.  (4.1) 

From  (4.1)  we  see  that  the  solution  can  be  obtained  along  each  ti  =  constant  plane  independent  of 
other  Ti.  In  particular,  this  implies  the  solutions  obtained  in  the  ii  <  0  half  space  can  be  obtained 
independently  of  the  solutions  in  the  -q  >  0  half;  by  the  symmetry  of  the  problem,  no  new 
solutions  will  be  found  in  the  ti  <  0  half  space  that  are  not  found  in  the  upper.  We  therefore 
restrict  our  discussion  to  the  half  space  t]  >  0  (y  >  x).  The  solution  in  the  plane  -r)  =  0  is 
discussed  in  Lemma  4.4  . 

In  the  -q  >  0  half  space,  the  solution  in  each  r\  =  constant  plane  is  just  that  for  a  generalized 
Riemann  problem  in  one  dimension. 

Lemma  4.2:  Let  u{t,  5,  tiq)  denote  the  generalized  Riemann  problem  solution  in  the 
Ti  =  Tio  >  0  plane.  Then  the  solution  «(f,  5,  tij)  to  the  generalized  Rionann  problem  in  the  nj 
plane  (tij  >  0)  is  just  a  sii  ilarity  solution  of  u(t,  ^,  tjq)  given  by 

m(c  t,ci  +  6(c),  Til)  =  «(f,  5,  Tic)  .         _  (4.2) 

where  c  =  —  . 

Proof:  If  vo(^)  is  the  initial  data  for  the  generalized  Riemann  problem  on  t]q,  then 
**'c(t  ■  c  5  +  i(c))  =  Vo(5)  is  the  initial  data  for  the  generalized  Riemann  problem  on  t]^  where 
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^1 


c  =  —  and  b(c)  is  given  by  the  wedge  pattern  of  the  two  dimensional  Riemann  problem. 

Further,  the  function  u(t  -  cr,  {  -  c  ?  +  fe,  tij  -  c  tiq)  -  u(t,  ^,  tiq)  is  a  weak  solution  of  the 
generalized  Riemann  problem  on  •i]^  having  initial  data  wq  .  n 

Thus  the  solution  to  the  2-dimensional  Riemann  problem  (4.1)  ir  the  tti  >  0  half  space  is 
found  from  a  mapping  M:  R^  -  R^  which  is  a  similarity  transformation  ui  r,  ?,  ti  and  therefore 
continuous  onto  the  domain  <  >  0,  5i  ^  >  0  • 

Theorem  4.3:  The  unique  (in  the  sense  of  Knizkov)  solution  in  the  plane  t]  >  0  (y  >  x),  to 
(1.1)  with  initial  data  that  is  piecewise  constant  on  a  finite  number  of  wedges  focused  on  a  single 
point  in  the  plane,  with/i  -  /j  -  /,  /  €  C^  :  i?  -  ^,  /  having  at  most  one  inflection  point,  is 

a)  piecewise  smooth. 

b)  composed  of  non-linear  waves  and  constant  states  which  are  the  images  under  a  similarity 
map  M  defined  by  Lemma  4.2  of  the  rarefaction  and  shock  waves  and  constant  states  of  a 
generalized  Riemann  problem  in  one  dimension. 

c)  composed  of  curves  of  irregular  points  corresponding  to  images  under  the  map  M  of  the 
irregular  points  in  the  1 -dimensional  Riemann  problem. 

Proof:  Follows  immediately  from  Lanma  4.2  and  Theorem  3.1  .  n 

Lemma  4.4: 

a)  If  the  half  line  -n  =  0,  5  <  0  is  not  a  line  of  discontinuity  of  the  initial  data  (ie:  is  not  a 
wedge  line)  the  waves  incident  upon  the  corresponding  half  plane  r  >  0,  ^  <  0,  t]  =  0  are 
continuous  across  the  half  plane. 

b)  K  the  half  line  ti  =  0,  ^  <  0  is  a  line  of  discontinuity  (constant  juir^j)  of  the  initial  data,  it 
remains  a  half  plane  t  >  0,^<0,y\  =  0  of  (in  general  variable)  jump  discontinuity  in  the 
solution. 

Proof: 

a)  Follows  from  taking  the  limit  ti  -  0  of  the  solution  in  either  half  plane. 

b)  Follows  from  taking  the  limit  of  solutions  for  which  the  wedge  line  in  the  initial  data  is 
displaced  slightly  from  the  ti  =  0  axis.  □ 

Similar  statements  hold  for  the  half  line  ti  =  0,  ^  >  0. 
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Although  the  map  M  provides  a  means  of  constructing  the  solution  to  the  2-dimensional 
Ri^nann  problem,  it  is  not  a  convenient  method  of  doing  so.  In  particular  it  is  a  method  that  will 
not  generalize  to  cases  f  *  g  where  no  appeal  to  a  one  dimensional  analysis  can  be  made.  We 
propose  that  the  correct  method  of  dealing  with  the  general  solution  to  the  2-dimensional  Riemann 
problem  is  to  identify  the  general  two  dimensional  non-linear  waves.  In  [4]  we  identify  these 
waves  for  the  case  /  ■  g  and  use  the  construction  method  initiated  by  Guckenhdmer  [1]  and 
Wagner  [6]  to  piece  these  waves  together  into  entropy  obeying  solutions. 

5.  Conduding  Remarks 

The  condition /i  ■  /2  ■  /,  f  (.  C^  :R  -R,  f  having  at  most  one  inflection  point  has  been 
shown  to  be  sufficient  to  guarantee  that  the  solution  to  the  2-dimensional  Riemann  problem  is 
piecewise  smooth.  A  construction  has  been  presented  detailing  a  presumed  mechanism  for  loss  of 
piecwise  smoothness  if  /  has  three  or  more  inflection  points. 

The  Riemann  problem  in  one  dimension  has  been  generalized  to  include  a  finite  number  of 
discontinuities  and  smooth  'rarefaction*  variation  in  the  initial  data.  Sufficiency  conditions  based 
upon  the  number  of  inflection  points  in  the  convection  function  /  have  been  obtained  wtadi 
guarantee  piecwise  smoothness  of  the  solution  to  the  generalized  problan. 
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Figure  CapdoDS 

Fig.  2.1  a)  The  upper  convex  envelope  drawn  between  two  points  u,  and  u^  for  an  illustrative 
function  /.  b)  The  Rieraann  solution  in  the  t,x  plane  based  on  a).  Dark  lines 
represent  shoci:  waves,  light  lines  are  rarefaction  waves. 

Fig.  2.2  A  sketch  of  the  solution  (lower  figure)  for  a  particular  case  of  initial  data  (upper 
figure)  for  a  generalized  Riemann  problem  with/  having  no  inflection  point.  The  dark 
line  is  the  shock  wave.  Important  rarefaction  lines  have  been  labelled  by  their  speeds. 

Fig.  2.3  A  sketch  of  the  solution  (lower  figure)  for  a  particular  case  of  initial  data  (upper 
figure)  for  a  generalized  Riemann  problem  with/  having  a  single  inflection  point.  The 
dark  line  is  the  shock  wave.  Important  rarefaction  lines  have  been  labelled  by  their 
speeds. 

Fig.  2.3  Formation  of  a  single  'bubble'  in  a  solution.  The  dotted  lines  in  the  middle  diagrams 
on  the  right  emphasise  the  non  -  linearity  of  the  convex  envelope. 

Fig.  3.1  Two  possibilities  encountered  under  case  P3  described  in  the  text.  Dark  lines  are  shock 
waves,  light  lines  are  rarefaction  waves. 
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